1 Astronomy 

(DOI: will be inserted by hand later) 
© ESO 2008 



Astrophysics 



o 
o 

(N 
G 



NLTE determination of the sodium abundance in a homogeneous 
sample of extremely metal-poor stars 

S.M. Andrievsky 1 ' 2 , M. Spite 1 , S.A. Korotin 2 , F. Spite 1 , P. Bonifacio 1 - 3 , R. Cayrel 1 , 

V. Hill 1 , and P. Francois 1 

1 GEPI, CNRS UMR 8111, Observatoire de Paris-Meudon, F-92125 Meudon Cedex, France 

2 Department of Astronomy and Astronomical Observatory, Odessa National University, Shevchenko Park, 65014 Odessa, Ukraine 

3 CIFIST Marie Curie Excellence Team 

ABSTRACT 
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Context. Abundance ratios in extremely metal-poor (EMP) stars are a good indication of the chemical composition of the gas in the earliest 
phases of the Galaxy evolution. It had been found from an LTE analysis that at low metallicity, and in contrast with most of the other elements, 
the scatter of [Na/Fe] versus [Fe/H] was surprisingly large and that, in giants, [Na/Fe] decreased with metallicity. 
^» . Aims. Since it is well known that the formation of sodium lines is very sensitive to non-LTE effects, to firmly establish the behaviour of the 
' sodium abundance in the early Galaxy, we have used high quality observations of a sample of EMP stars obtained with UVES at the VLT, and 
^\ , we have taken into account the non-LTE line formation of sodium. 

t-H ■ Methods. The profiles of the two resonant sodium D lines (only these sodium lines are detectable in the spectra of EMP stars) have been 
' computed in a sample of 54 EMP giants and turn-off stars (33 of them with [Fe/H] < -3.0) with a modified version of the code MULTI, and 
, compared to the observed spectra. 

Results. With these new determinations in the range -4 < [Fe/H] < -2.5, both [Na/Fe] and [Na/Mg] are almost constant with a Jow scatter. 
In the turn-off stars and "unmixed" giants (located in the low RGB): [Na/Fe] = -0.21 ±0.13 or [Na/Mg] = -0.45 ±0.16. These values 
C"| are in good agreement with the recent determinations of [Na/Fe] and [Na/Mg] in nearby metal-poor stars. Moreover we confirm that all 
Oh the sodium-rich stars are "mixed" stars (i.e., giants located after the bump, which have undergone an extra mixing). None of the turn-off 
^ , stars is sodium-rich. As a consequence it is probable that the sodium enhancement observed in some mixed giants is the result of a deep mixing. 
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' 1. Introduction Such abundance ratios as [a/Fe], [Na, Al/Fe], and others can 

h3 . be derived from high-resolution spectra and then compared to 

.«? , During the last decade the interest in the extremely metal-poor the SNe n yMd prediction In particu i ar m Cayrel et al. (2004) 

(EMP) galactic stars has significantly increased. These stars, in and in Bonifacio et al (2006X a homog eneous sample of EMP 

fact, are reservoirs of the primordial Galactic material that was ^ has been studied . 35 giants and 18 stars (39 of 

only slightly polluted by nuclides produced in SNe II (by the them wkh [Fe/H] < _ 3) It has been shown ^ fa these EMp 

first generation of massive stars). Thus, elemental abundances ^ ^ scatter of ^ abundanc e ratios at metanicity be i ow 

provided by metal-poor old stars, and especially abundance ra- _ 3 q dex is generally very small and comparable to the mea- 

tios, are of the highest importance for the testing of the SN surement errors But ±eK are some exce ptions, and one of the 

theoretical models, their yields, and finally for constraining the mofjt striking casefj ifj [Na/Fe]> which varies from star t0 star 

initial mass function slope in the early Galaxy and understand- by a factor dose tQ ten In these EMp stars> only the resonance 

ing its formation history and early evolution. lines of Na (D1 and m) are visibl ^ and these lines are known 

With the advent of large aperture telescopes, the precise el- to be very sensitive to NLTE effects (e.g., Mashonkina et al. 

emental abundances in these faint objects became available. 1993; Baumiiller et al. 1998). 



Send offprint requests to: M. Spite 
e-mail: Monique.Spite@obspm.fr 



In Cayrel et al. (2004) and Spite et al. (2005) the sodium 
abundance has been calculated under the LTE hypothesis. 
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Fig. 1. [Na/Fe] vs. [Fe/H] for an homogeneous sample of EMP 
stars (Cayrel et al. 2004; Bonifacio et al. 2006): dwarfs (black 
dots), low RGB stars (grey dots), and upper RGB stars (open 
circles). The abundance of Na is deduced from the D lines 
through an LTE analysis. As a mean, the ratio [Na/Fe] is lower 
in turn-off stars (black dots) than in RGB stars and is also 
lower in lower RGB stars (grey dots) than in upper RGB stars 
(open circles). This graph suggests a variation of the computed 
sodium abundance with the gravity, which could be the conse- 
quence of an NLTE effect. 



Later, as a very first approximation, a uniform correction LTE- 
NLTE of -0.5 dex, extrapolated in the table of Baumiiller 
et al. (1998) for log g — 2, was applied to all the giants. 
The scatter of the relation [Na/Fe] versus [Fe/H] in giants 
was then very large, and [Na/Fe] was decreasing with the 
metallicity in contrast to the turn-off stars (Spite et al., 2005) 
where the ratio [Na/Fe] was almost constant in the interval 
-3.5 < [Fe/H] < -2.5. 

In Fig. [T] we show the relation [Na/Fe] vs. [Fe/H] for our 
sample of dwarfs and giants without any NLTE correction. The 
lower RGB stars (below the bump, see Spite et al. 2006) and 
the upper RGB stars (above the bump) have been plotted with 
different symbols. In the mean, the dwarfs have a lower value 
of [Na/Fe] than the giants, and the lower RGB stars a lower 
value of [Na/Fe] than the upper RGB stars. This variation of 
[Na/Fe] with the gravity of the star suggests that, as expected, 
non-LTE effects are very pronounced on the Na D lines, and 
that as a consequence it is important to carefully take them 
into account. In this paper we present new determinations of 
the sodium abundance in the sample of EMP stars of Cayrel et 
al. (2004) and Bonifacio et al. (2006) through a careful NLTE 
analysis, and we discuss the new trends obtained. 



2. The star sample 

The spectra used here have been presented in detail in Cayrel 
et al. (2004) and Bonifacio et al. (2006). The observations were 
performed with the ESO VLT and its high resolution spectro- 
graph UVES (Dekker et al. 2000). The spectra have a resolv- 



ing power of 43000 in the region of the D lines and a typi- 
cal S/N ratio per pixel of ~ 200; they have been reduced us- 
ing the UVES context (Ballester et al. 2000) within MIDAS. 
Stellar parameters (effective temperature, surface gravity, mi- 
croturbulent velocity, and metallicity) for each star have been 
taken from Cayrel et al. (2004) and Bonifacio et al. (2006). 

3. Method of analysis 

To derive the sodium abundance in the program stars, we per- 
form an NLTE analysis of the two resonant sodium Dl and D2 
lines (only these lines are seen in the spectra of the EMP stars) 
using a modified version of the MULTI code (Carlsson 1986). 
These modifications are described in Korotin et al. (1999ab). 
To more adequately calculate the continuous opacity taking 
into account the absorption produced by the great number of 
spectral lines, the additional opacity sources from ATLAS9 
(Kurucz 1992) are included. Simultaneous solution of the ra- 
diative transfer and statistical equilibrium equations are ob- 
tained using the approximation of a complete frequency redis- 
tribution for all the sodium lines. Atmosphere models of appro- 
priate metallicity are interpolated in Kurucz's grid of stellar at- 
mosphere models (microturbulent velocity 2 km s , a — 1 .25). 
This modified code has already been used for NLTE abundance 
determination of carbon, oxygen, and sodium in stars of differ- 
ent types (see, e.g., Korotin et al. 1999ab; Korotin & Mishenina 
1999; Andrievsky et al. 2001; Andrievsky et al. 2002, and ref- 
erences therein). 

We employ the modified model of the sodium atom (first 
considered by Sakhibullin 1987) that consists of 27 energy lev- 
els of Na I atoms in addition to the ground level of the Na II ion. 
The fine splitting has been taken into account only for the 3p 
level. This enables one to calculate the sodium doublet more 
precisely. The radiative transitions between the first 20 levels 
of Na I and the ground level of Na II are considered, while 
transitions between the other levels are used for the particle 
number conservation. Linearization includes 46 b - b and 20 
b — f transitions. Radiative rates for 34 transitions are fixed. 
Photoionization rates were taken from TOPBASE. Collisions 
with electrons, as well as with hydrogen atoms, have been in- 
cluded. 

Figure |2] shows the behaviour of the departure coefficients 
of Na I for dwarfs and giants of different metallicities. It is 
clearly seen that within the region of the line-formation in the 
atmosphere, the lower levels of the Na I atom are overpopu- 
lated, leading to a larger opacity in Dl,2 lines compared to 
the pure LTE case, and therefore to larger equivalent widths 
of these lines. This means that NLTE abundance corrections 
will be negative. 

To find the relative-to-solar sodium abundance, we com- 
puted NLTE sodium abundance in the Sun. This was done 
with the following lines: 4496.05, 4982.81, 5148.84, 5682.63, 
5688.20, 5889.95, 5895.92, 6154.22, 6160.75, 8183.25, 
8194.82, and 12679.14 A. Their profiles were extracted from 
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Fig. 2. Departure coefficients in atmospheres of dwarfs and gi- 
ants of the different metallicities. 




5890.0 5890.5 A (A) 



5896.5 I (A) 



Fig. 3. The observed spectra of the Dl and D2 lines (dots) 
in BS 16477-003 are compared to the theoretical profiles 
computed for the best-fit Na abundance +0.10 and - 0.10 
(solid lines). 



the Kurucz et al. (1984) solar flux spectrum. Similarly to 
the stellar spectra analysis, to derive the NLTE sodium abun- 
dance in the Sun, we used the solar atmosphere model from 
the Kurucz's grid and the solar microturbulent velocity rec- 
ommended by Maltby et al. (1986). Our NLTE solar sodium 
abundance is e(Na) = 6.25 + 0.04, in agreement with the 
value found by Baumuller et al. ( fl998"l >: e(Na) = 6.30 + 0.03, 
and Mashonkina et al. (EOOOl l: e(Na) = 6.20. The value 
e(Na) Q = 6.25 will be used later, as a reference to determine 
[Na/H]. 



4. Results 

New sodium abundances have been deduced from NLTE pro- 
files of the D lines for our sample of EMP stars. Typical un- 
certainty in the derived sodium abundance is about +0.05 (see 




5890.0 5890.5 A. (A) 



5895.5 5896.0 5896.5 A (A) 



Fig. 4. Best fit between observed (dots) and calculated 
(solid line) profiles of the Na Dl (5889.95 A) and D2 
(5895. 92A) lines, for some program stars. 



Fig. |3]l. Figure |4] shows the best fitting between observed and 
calculated profiles for some EMP giants and dwarfs. In Table 
1 the parameters of the models (Teff, log g and v t ) which have 
been used and the new sodium abundance in logarithm e(Na) 
(for e(H) = 12) are given for each star. The stellar ratios [Na/H] 
have been computed assuming e(Na) = 6.25 (see Sect. [3). 

Three stars in our sample are very peculiar carbon- 
rich stars: CS 22892-52 (Sneden et al. 1996, 2000, 2003), 
CS 22949-37 (Norris et al. 2001; Depagne et al. 2002), and 
CS 29528-41 (Sivarani et al. 2006). The models we used for 
the NLTE computations do not take the carbon enhancement in 
these stars into account and thus the sodium abundance com- 
puted here must be considered as very preliminary. These stars 
are not taken into account in the discussion. 

In Fig.[5]we present an estimate of the difference in sodium 
abundance obtained when LTE or NLTE profiles of the D lines 
are used along the Red Giant Branch of the HR diagram. The 
computations have been done for a set of temperatures between 
6500 K (turn-off stars) to 4500 K (upper RGB stars) and metal- 
licities from -2.5 to -4.0 (Table 2). For each temperature the 
gravity has been deduced from the isochrones of Kim et al. 
(2002) for 14 Gyr, [a/Fe] = +0.6, [Fe/H] = -3.7and - 2.7 
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Table 1. Model parameters and NLTE sodium abundance in our sample of stars. The three stars marked with an asterisk are 
carbon-rich stars. An "m" in the last column denotes the mixed giants following Spite et al. (2005). 
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Table 2. Abundance corrections (e(NLTE) - e(LTE)) for the 
Dl and D2 lines. 
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Fig. 5. Individual correction applying to an LTE analysis of the 
Na Dl and D2 lines along the RGB, for metal-poor stars with 
-4.0 < [Fe/H] < -2.5 and the hypothesis that [Na/Fe] = . 



with the classical formula: 



log L/Lq = log 



+ 4 log T eff /T eff>0 - log g/g 



For this set of models, the LTE profiles have been com- 
puted with the hypothesis that [Na/H] = [Fe/H] (i.e., [Na/Fe] = 
0), and then fitted by NLTE profiles that enabled the derivation 
of the appropriate NLTE corrections. A first cause of uncer- 
tainty is that the shapes of LTE and NLTE profiles are slightly 
different, but the correction mainly depends on the equivalent 
width of the sodium lines, as can be seen in Fig. [6] As a con- 
sequence, the corrections given in Table [2] can be used only if 
the equivalent width of the lines is similar to that given in the 
table. If there is a strong anomaly of the sodium abundance 
in the star, it is better to deduce the correction from the Fig. 
[6] (or to directly compute the non-LTE profile of the line as 
it has been done for the stars of our sample). Takeda et al. 
(2003) have made a first attempt to compute the non-LTE cor- 
rection of the sodium abundance in a large metallicity range 
(-4.0 < [Fe/H] < +0.4). We have checked that in our domain 
of metallicity ([Fe/H] < -2.5) there is a rather good agreement 
with their computations. 
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5. Discussion 

5.1. Mixing effects in giant stars 

In Fig.|7]we present the NLTE values of [Na/Fe] for our sample 
of stars as a function of [Fe/H] (the carbon-rich stars have been 
omitted). The symbols are the same as in FigQ] Dwarfs and gi- 
ants now lie in the same locus, and the trend of the [Na/Fe] ra- 
tio, in the range -4.0 < [Fe/H] < -2.5, is flat. The mean value 
of [Na/Fe] is close to -0.2. 

A few stars appear to be "Na-rich" with [Na/Fe] > +0.1. 
All these stars belong, following Spite et al. (2006), to the class 
of the "mixed stars". They present the characteristics of a mix- 
ing with the H-burning shell, low abundance of carbon corre- 
lated with a high abundance of nitrogen, and a low value of 
the 12 C/ I3 C ratio close to the equilibrium value. However, all 
the mixed stars are not Na-rich. Let us remark that all the stars 
found here with [Na/Fe] > +0.1 were already suspected from 
LTE computations, to be Na-rich. All of them, but BD 17:3248, 
have also been found to be aluminum rich, and from their po- 
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Fig. 6. Individual corrections plotted as a function of the equiv- 
alent widths of the lines for the different models given in Table 
[2] The solid line represents the corrections for models with 
[Fe/H]=-2.5, the dashed lines models with [Fe/H]=-3.0, the 
dotted-dashed lines models with [Fe/H]=-3.5, and the dotted 
lines models with [Fe/H]=-4.0. For a same equivalent width of 
the D lines the correction little depends on the parameters of 
the models 

sition in the HR diagram they could be AGB stars (see also 
Johnson d2002l and Bond ([19801 ) for BD 17:3248). None of the 
turn-off stars or "unmixed" giants are sodium rich. As a con- 
sequence the sodium rich stars are suspected to have suffered 
a mixing between the atmosphere and the H-burning shell, a 
mixing deep enough to bring the products of the Ne-Na cy- 
cle to the surface. When only EMP dwarfs and unmixed giants 
are taken into account, we find that [Na/Fe] = -0.21 +0.13 
or [Na/Mg] = -0.45 + 0.16; the star-to-star variation becomes 
small and comparable to what is obtained for the other elements 
(see Cayrel et al. 120041 

5.2. Comparison with the nearby metal-poor stars, 
and the prediction of the models of galactic 
chemical evolution 

Recently Gehren et al. (2006 ) have computed the sodium abun- 
dance in a sample of 55 nearby metal-poor stars, most of them 
belonging to the thin or thick Galactic disc and some of them to 
the halo. They have very carefully taken into account the non- 




[Fe/H] 



Fig. 7. [Na/Fe] vs. [Fe/H] for the sample of EMP stars. The 
symbols are the same as in Fig.Q] The sodium abundance has 
been computed in NLTE with a modified version of the pro- 
gram of Carlsson. The trend of the [Na/Fe] ratio is flat in the 
region -4.0 < [Fe/H] < -2.5, and there is a good agreement 
between the abundances of the giants and turn-off stars. All the 
stars with [Na/Fe] > 0.1 have been found to be extra-mixed, in 
Spite et al. d20l)5l 



LTE effects. In Fig.[8]we compare the ratio [Na/Fe] in our EMP 
stars and in the sample of Gehren et al. ([2006 ); unfortunately, 
we have no star in common to check the influence of system- 
atic effects, but they have one star in our domain of metallicity 
(G64-12 with [Fe/H] =-3.1), and this star falls exactly in the 
middle of our sample in Fig. [8] 

Sodium has only one (stable) isotope, 23 Na, which is 
mainly produced in the carbon burning process operating in- 
side massive stars. As noted by Woosley & Weaver (1995) the 
resulting amount of sodium nuclei in the SN ejecta is sensi- 
tive to the excess of neutrons, therefore its production is some- 
what metal dependent, but in the early Galaxy (very low metal 
content) this element production should be the same as for the 
primary elements. There were several attempts to incorporate 
the observed sodium abundance in the disc and halo stars into 
Galactic evolutionary models. One can mention here the works 
of Timmes et al. (1995), Samland (1998), and Goswami & 
Prantzos (2000). All the above-mentioned models are based on 
the Woosley & Weaver (1995) yields from massive stars, but 
use different initial-mass functions, star formation rates, and 
other model specifications, which can explain the slight differ- 
ences in the predictions. 

In Fig.[8]our observations and those of Gehren et al. (2006) 
are compared to the predictions of Timmes et al. ( I19951 ). These 
authors predict a local minimum around [Fe/H] = -1.5 with 
subsequent increasing of [Na/Fe] as the metallicity decreases 
to [Fe/H] = -3. The predictions for two different values of the 
iron yields from massive stars (which depend on the mass cut in 
the massive supernovae) are represented in the figure. Between 
models T(A) and T(B), the mass of iron ejected by SN II dif- 
fers by a factor of two. At a metallicity higher than [Fe/H] = 
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Fig. 8. [Na/Fe] vs. [Fe/H] for the sample of EMP stars (Cayrel 
et al. 120041 Spite et al. 120051 ), dwarfs (black dots), unmixed 
giants (grey dots), and the nearby metal-poor stars studied by 
Gehren et al. (2006) (open triangles). The dashed lines repre- 
sent the predictions of two models of Timmes et al. ( 1995 ) that 
differ by the quantity of iron ejected by massive stars and thus 
by the position of the mass cut. 



-3, there is a rather good agreement with the predictions of 
the model T(A), which supposes the smallest value of the iron 
yields from massive stars, and such an agreement is generally 
found for the other elements (see Cayrel et al. 120041 1. 

In contrast with [Na/Fe], the ratio [Na/Mg] is independent 
of the mass cut. In Fig. [9] we present [Na/Mg] as a function of 
[Fe/H] for our sample of EMP stars and the sample of Gehren 
et al. d20061 l. (We neglected the NLTE effects on the magne- 
sium lines that are not supposed to be very important). The 
agreement is rather good for [Fe/H] > -3. 

Timmes et al. ( 1 19951 1 do not predict the behaviour of 
[Na/Fe] or [Na/Mg] for [Fe/H] < -3. From our measurements, 
[Na/Fe] and [Na/Mg] become flat: the relative abundances are 
independent of the metallicity. However, a slight increase of 
[Na/Mg] when the metallicity decreases (from [Na/Mg] = -0.5 
at [Fe/H]=-3. to [Na/Mg] = -0.3 at [Fe/H]=-4.) cannot be ex- 
cluded. 

Samland (1998) has considered the chemical evolution of 
many elements and found that the expected distribution of 
[Na/Fe] vs. [Fe/H] in the extreme metallicity region from 
[Fe/H] = -4 to -2.5 is almost independent of [Fe/H], and 
the mean [Na/Fe] value equals approximately -0.2. This was 
not supported by observations available at that time, but is 
in good agreement with our observations. Qualitatively the 
same conclusion was reached by Goswami & Prantzos (2000). 
According to their model, [Na/Fe] must be almost constant 
within this metallicity interval and equal approximately -0.5, 
but this is significantly lower than what is given by our calcu- 
lations. 

Recently Tsujimoto et al. (120021 1 have also investigated the 
chemical evolution of sodium and aluminum during the early 
epochs of the Galaxy. They have used a model that reproduces 




[Fe/H] 

Fig. 9. [Na/Mg] vs. [Fe/H] for the sample of EMP stars and 
the nearby metal-poor stars studied by Gehren et al. (2006). 
Symbols are the same as in Fig. [8] The dashed line repre- 
sents the prediction of the model of Timmes et al. (1995). 
(The [Na/Mg] ratio is independent of the mass cut at vari- 
ance with [Na/Fe] and the predictions are reduced to a sin- 
gle curve). It can be seen in the figure that in the interval 
-4 < [Fe/H] < -2.8 [Na/Mg] is constant. However, a slight 
increase of [Na/Mg] when the metallicity decreases (from 
[Na/Mg] = -0.5 at [Fe/H] =-3. to [Na/Mg] = -0.3 at [Fe/H]=- 
4.) cannot be excluded. 



the distribution of EMP stars in the [Na/Mg] or [Al/Mg] ver- 
sus [Mg/H] plane. The observed trends and the scatter of the 
observations are mainly based on the work of Mc William et 
al. ( 119951 ). This sample contains several mixed giants which 
explain the very large scatter of [Na/Mg] and [Al/Mg] they 
observed at low metallicity. Such a scatter, which is an im- 
portant point in the definition of the model of Tsujimoto et 
al. d2002 ), is not observed in the unmixed stars of our sam- 
ple (dwarfs and "unmixed" giants) and thus does not probably 
reflect a scatter in the pristine material. Let us note that Cohen 
et al. (120041 ). who studied a sample of turnoff stars in the in- 
terval -3.6 < [Fe/H] < -2.2, also found a very small scatter of 
[Al/Mg]. The precise determination of [Al/Mg] in our sample 
of stars from non-LTE synthetic profiles is in progress. 

6. Conclusion 

NLTE abundances of sodium have been derived for 
a sample of extremely metal-poor dwarfs and giants us- 
ing high-resolution spectra obtained with VLT. Both sub- 
samples of dwarfs and unmixed giants now have the 
same mean [Na/Fe] or [Na/Mg] values, within the standard 
errors: [Na/Fe] = -0.21 ±0.13 or [Na/Mg] = -0.45 + 0.16. 
However, a slight increase of [Na/Mg] at low metallicity, (from 
[Na/Mg] = -0.5 at [Fe/H]=-3. to [Na/Mg]= -0.3 [Fe/H]=^k) 
cannot be excluded. 

• In the domain -4 < [Fe/H] < -2.5, the giants of our sam- 
ple do not exhibit any systematic dependence of [Na/Fe] upon 
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[Fe/H], contrary to the preliminary findings of Cayrel et al. 
(2004). This effect was a consequence of a coarse treatment 
of the non-LTE effects: in Cayrel et al. (2004) the LTE abun- 
dances had been uniformely corrected by -0.5dex. In fact, the 
non-LTE correction strongly depends on the equivalent width 
of the sodium lines and thus is generally smaller at lower metal- 
licity. (When the lines are smaller they form deeper in the at- 
mosphere where, in particular, collisions are more important 
and thus non-LTE corrections are smaller.) 

• Among the giant stars, the "mixed" giants (stars that have 
undergone an extra-mixing between the atmosphere and the 
H-burning shell following Spite et al. 2005 ), often present an 
overabundance of Na. It is suggested that in these stars, mate- 
rial processed in Ne-Na cycle in the hot hydrogen burning shell 
has been brought to the surface. From their position in the HR 
diagram, some of these stars could be AGB stars. 

• Qualitatively, our observational result on sodium distri- 
bution in the metallicity region [Fe/H] from -4 to -2.5 agrees 
with predictions of models elaborated by Samland (1998) and 
Goswami & Prantzos (2000). Both models produce the al- 
most constant [Na/Fe] value in the region of extreme metal- 
poor stars. In particular, the Samland's model gives an average 
[Na/Fe] value that approximately equals -0.2, and this value 
appears to be in quantitavive agreement with our estimate of the 
mean relative-to-iron sodium abundance in the early Galaxy. 
The model of Timmes et al. ( 1995 ) also agrees with our obser- 
vations around [Fe/H] =-3. In qualitative agreement, the model 
of Goswami & Prantzos predicts a slightly too low relative 
sodium abundance for extremely metal-poor stars. 
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